arXiv: 1509.05047vl [astro-ph.GA] 16 Sep 2015 


Draft version September 21, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


ON THE IMF IN A TRIGGERED STAR FORMATION CONTEXTt 
Tingtao Zhou^’^, Chelsea X. Huang^’^, D.N.C. Matthias Gritschneder'*’®, & Herbert Lau'’ 


1 


Kavli Institute for Astronomy &; Astrophysics and School of Physics, Peking University, Beijing China, edmondztt@gmail.com 

^Department of Physics, MIT, USA 
^Department of Astrophysical Sciences, Princeton University, USA 
"^UCO/Lick Observatory, University of California, USA 
^Institute for Advanced Studies, Tsinghua University, Beijing, China. 

^University Observatory Munich, Germany and 
^Argelander Institute, University of Bonn, Germany 
Draft version September 21, 2015 


ABSTRACT 

The origin of the stellar initial mass function (IMF) is a fundamental issue in the theory of star 
formation. It is generally fitted with a composite power-law. Some clues on the progenitors can be 
found in dense starless cores which have a core mass function (CMF) with similar shape. In the low- 
mass end, these mass functions increase with mass, albeit the sample may be somewhat incomplete, 
whereas they decrease with mass in the high-mass end. There is an offset in the turn-over mass between 
the two mass distributions. The stellar mass for the IMF peak is lower than the corresponding core 
mass for the CMF peak in the Pipe nebula by about a factor of three. Smaller offsets are found between 
the IMF and the CMFs in other nebulae. We suggest that the offset is likely to be induced during a 
star-burst episode of global star formation which is triggered by the formation of a few O/B stars in the 
multi-phase media, which is naturally emerged by the onset of a thermal instability in the cloud-core 
formation process. We consider the scenario that the ignition of a few massive stars photoionizes the 
warm medium between the cores, increases the external pressure, reduces their Bonnor-Ebert mass 
and triggers the collapse of some previously stable cores. We quantitatively reproduce the IMF in the 
low-mass end with the assumption of additional rotational fragmentation. 

Subject headings: stars: formation-ISM: individual objects: Pipe Nebula - HII regions-ISM: clouds- 
ISM: structure - methods: analytical 


1 . INTRODUCTION 

Recent infrared measurements of dust extinction, as 
well as CO and NH 3 maps of the filaments in molecu¬ 
lar clouds reveal a population of embedded cores mostly 
confined by the global external pressure from the inter¬ 
core gas. These cores are closely associated with the pro¬ 
genitors of young stellar objects. In such filaments there 
are gravitationally bou nd cores (or pre-ste llar cores), and 
embedded protostars (jAndre et al.ll20I(i|l . But, in the 
Pipe nebula nearly all the cores are starless (iLada et al.l 
l2008f) . Only the most massive cores in Pipe are gravi¬ 
tationally bound and migh t collapse. The s hape of the 
CMF of these dense cores (ILada et al.l 1200811 appears to 
be qualitati vely similar to the broken power-la w of the 
stellar IMF (jKrouDall2002l: IWeidner et ani2nilh . An in- 
depth analysis of their structure and evolution may be 
useful for the construction of a star formation theory. 

One noticeable difference between the CMF and the 
IMF is an offset between the two distributions. The ra¬ 
tio of the characteristic mass of the IMF over the char¬ 
acteristic mass of the CMF is generally smaller than 
unit y. This ratio ranges fro m ~ I : 3 in the Pipe neb- 
ula (iRat hborne et al. 12006 1 to > I : 2 in the Orion 
( Nutter fc Ward-T hom pson 1200711 and Aquil a nebulae 
( Konvves et al.ll2010D 7 lRathborne et al.l (1200911 proposed 
that this offset is due to a direct one-to-one mapping from 
the cores to the stars with sufficient amount of mass loss 
during the star formation process. It is hard to explain 
in this scenario how the low-mass stars are produced by 
the originally stable low-mass cores. 


Many studies modeling the origin of IMF highlight 
factors such as the accretion rate of protostars, tur- 
bulent fragmentation and accretion of c ores (see e.g., 
Bonnell fc Bate! 120061 : IClark et'al] 120081 : lAnathnindil^ 
2 OI 1 I 1 . In the meantime, constraints such as the small 
age spread of stars in young clusters are inconsistent in 
these models. 

We propose a scenario to explain the transition from 
the CMF to the IMF as well as a synchronized star forma¬ 
tion, triggered by a few fi rst formed 0/B stars in the neb- 
ula. Many authors ('e.g.. IStromgrenl[l939l: ISpitzedllOTSU 
studied the effects of UV radiation from massive stars 
onto the surrounding regions. The UV radiation gen¬ 
erates a hot ionized region (Stromgren sphere), increas¬ 
ing the ambient pressure. Consequently, an isothermal 
shock is driven through the nebula. Influenced by the 
Stromgren sphere, denser sub-structures are rapidly en¬ 
hanced and most of the pre -existing pressure confine d 
cores are compressed (e.g., iGritschneder et all 1200911 . 
The ionization timescale is generally short compared to 
the hydrodynamical timescales, such that the increase in 
the background pressure and temperature can be consid¬ 
ered instant in most cases. 

In this work, we develop a quantitative understanding 
of the consequences of this sudden change in the ambi¬ 
ent pressure and temperature in the context of the tran¬ 
sition from the CMF to the IMF. A direct result is the 
reduction of t he critical Bonnor-Ebert mass (|Ebertill955l : 
IBonnod 1195611 . making previous unbound cores gravita¬ 
tionally unstable. Therefore, the whole region experi¬ 
ences a rapid star burst, synchronizing star formation. 
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The shifted Bonnor-Ebert mass also naturally leads to 
the transition from the CMF distribution to the IMF 
distribution. 

We base our calculations on a Pipe-like cloud, with 
a starless CMF as the starting point of the evolution. 
We investigated the formation of a starless CMF in our 
previous work (|Huang et al.l l20ll . hereafter HI3). In 
H13, we suggest that, prior to the onset of global star 
formation, the cores and the inter-core gas are two sepa¬ 
rate phases in pressure eq uilibrium, which poss ibly result 
from thermal instabili ty (iLin fc Murra^l200flfl . The dy- 
namics of coagulatio n (jMnrrav fc Linlll99lH) and ablation 
(jMurrav et al.l [199^ dominate the evolution of the sys¬ 
tem and the star formation timescale is prolonged by the 
turbulent pressure or magn etic pressure inside the cloud 
(jLazarian fc Vishniaclll99^ . 

In this paper, we consider the consequences of the ig¬ 
nition of the first massive stars in a starless nebula. In 
§ 2 , we suggest that the photoionization of the inter-cloud 
gas leads to a decrease in the critical Bonnor-Ebert mass 
and triggers a global burst of star formation. We use 
this reduced critical mass and the starless CMF from the 
Pipe to obtain the transition to a stellar IMF in §3. We 
assume different star formation statistics, showing that 
uncertainties such as binary mass ratio distribution do 
not affect the IMF shape in a certain parameter range. 
The discussion and conclusions are presented in §4. 

2 . TRIGGERED STAR FORMATION 

After the formation of the first massive stars, their 
stellar luminosity has a significant influence on the sub¬ 
sequent evolution of the surrounding media. In the pres¬ 
ence of a strong UV flux from nearby 0/B stars, the 
diffuse neutral medium, originally at 100 K, hereafter 
referred to as the warm medium, would be ionized and 
heated up, while the dense molecular region, originally 
at 10 K, hereafter referred to as the cold medium, would 
still be mostly shielded. The resulting increase in the ex¬ 
ternal pressure leads to a reduction in the Bonnor-Ebert 
mass of the cold cores and a decrease in the star forma¬ 
tion timescale. 


(last described by iFerland et al.l 1199811 . The effective 
temperature of the central star is taken to be 34,700 K 
with a surface luminosity of log(L/L 0 ) = 5.59. This is 
a typical value for a main sequence star with a mass of 
4 OM 0 , produ ced with the C ambridge stellar evolution 
code STARS (jEggletonl 1197111 . Assuming a blackbody, 
the star emits 31% of its power in hydrogen ionization 
photons, corresponding to ( 3 o ,49 1 . 

We assume a solar metallicity for the gas and a stan¬ 
dard dust composition. We use a constant hydrogen den¬ 
sity of riH = 774 cm“^ to simulate the warm medium 
in Pipe nebula, to be c onsistent with the values in 
iGritschneder fc Lliil (|2012[1 . 

The results show that there is a sharp temperature 
dropoff from several thousand degrees to several hundred 
degrees at a partic ular depth of the cloud. The solid 
line in Figure 1(a) shows the temperature of the warm 
medium versus the distance from the star. This indicates 
the Stromgren radius is 0.7 pc and the temperature inside 
the ionized sphere is around 9000 K. Given the typical 
size of a nebula, which is usually a few pc, only a few 
massive stars are needed to ionize most of the inter-core 
media, leading to a much higher exterior pressure. 


2.2. Heating of the cold cores 

The ionizing UV flux can also penetrate into dense 
cores. We first consider the heating from the radiation 
which balances the cooling from recombination. 

We denote i?ion to be the distance from an 0/B star 
at which a typical cold core is heated by the external 
radiation to 100 K. By neglecting the loss of ionizing 
photons from the star through the warm medium, we 
estimate i?ion from: 
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2.1. The ionization of the warm medium 

The ionizing photons from the native stars create 
Stromgren spheres in the surrounding warm medium. 
For a typical 0/B star with an effective UV photon 
{Ery > 13.6eV) emission rate of Qo 104 ®'^ 49 g-i^ 
Stromgren radius Rs is determined by 

47r 

Qo ~ -R| OtB (1) 

1 o in — 13 'T^O.SSS—0.034 log 74 3 —1 * .i-i-. 

where as = 2.6 x 10 1a ___ cm^s is the 

recombination coefficient (|Drainell2011ll . We define r 4 as 
the temperature of the warm medium in units of 104 K, 
Qo ,49 as the number of UV photons in units of 1049s“4^ 
and nwarm ,2 as the warm medium hydrogen number den- 

1—2 

sity in units of 10^ cm“^. Within i ?5 ~ 3 pc Qg 49 n^a™ 2 
of the ionizing sources, the warm medium is nearly fully 
ionized, with Twarm ^ 104K. This new temperature cor¬ 
responds to an increase by a factor of 10 ^ 100 . 

More realistic calculations of equilibrium temperature 
have been performed with version 08.00 of GLOUDY 


with as here is the same as in section nn 
More realistically, we compute the ionization in the 
cold core with GLOUDY, using the same setup as in 
§ 2 . 1 , but this time with an initial hydrogen density 
nil = 774 cm“^ in the range of 0 ^ 0.5 pc, while 
riH = 7300 cm“9 from 0.5 pc to faraway, corresponding 
to the cold medium placed 0.5 pc from the central star, 
and the warm medium filling the space in between. Then 
instead of keeping constant density, we keep the two me¬ 
dia in pressure equilibrium, as the central star illumi¬ 
nates its surroundings. In this case, the density of the 
cold medium is enhanced rapidly, due to the sudden in¬ 
crease of warm medium temperature and pressure. After 
the attenuation of the warm medium, the typical photon 
penetration depth inside the cold medium is only about 
0.02 pc. This is much smaller than the typical core radius 
Tc ~ 0.1 pc. Therefore, only cores very near the massive 
star or very small cores are subjected to this effect. The 
temperature inside the cores further away is limited by 
the cooling process to be around 15-20 K, which is two 
times higher than before the ignition of the UV flux (see 
Figure 1(b) for depth > 0.5pc). 
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In addition, we also investigate the classical evapo¬ 
ration of cores due to the high temperature environ¬ 
ment. Using the formula for the e vaporation rate of 
clouds embedded in a 10^ K gas (iGraham fc Langed 

ll973l : [Cowie fc McKeell977HMcKee fc Cowielll977ll . with 

the limitations that the background gas is either fully 
ionized or neutral, 

16tT^KT c fl.3 X 10^^T^/^rc,pcgS“^ , k = Kn 
25k ^ \2.75 X lO-^rs/^re pcgs-i, k = 

( 4 ) 

Here, T is the environment temperature, k is the Boltz¬ 
mann’s constant, Tc^pc is the radius of the core in units 
of pc, Kn is the neutral conductivity, Kc is the classi¬ 
cal conductivity for a fully ionized gas Q- For a typical 
core with mass me, the evaporation timescale therefore 
is Tevap.n = TOc(to)“^ = 8 X 10® r^^pc Myr in the neutral 
case and Tevap.i = 4 x 10®r^ peMyr in the fully ionized 
case. Based on these calculations, the classical evapora¬ 
tion only affects cores with radius as small as 0.01 pc. 

We conclude from the above calculation that the evap¬ 
oration effect is negligible. This is understandable as 
both the recombination timescale and cooling timescale 
inside the cold medium are much shorter than these 
timescales inside the warm gas due to their density con¬ 
trast. Therefore, the ionization fraction and temperature 
of the cold medium can be maintained at low levels. We 
do not include the change of the IMF due to the evapo¬ 
ration in the calculation below. 

2.3. Change in Bonnor-Ebert mass and star formation 

timescale 

The star f ormation times cale in general can be as long 
as 100 Myr (|Ostrikei] 1201 111 . However, the collapse pro¬ 
cess is speeded up considerably by th e compression of the 
dense cores due to the UV feedback (|Gritschneder et id] 
1200911 . More importantly, the enhancement of the exter¬ 
nal pressure reduces the Bonnor-Ebert mass and induces 
the formation of low-mass stars. 

The Bonnor-Ebert mass of a cold core can be expressed 
as mBE oc wh ere Payt, is the ex ternal pressure 

(see also equation(4) in lLada et al.ll200^ . Prior to the 
formation of a massive star, the temperature of the warm 
medium is Text ~ 100 K, and the cores’ temperature is 
Tint 10 K. In a pressure equilibrium, the density con¬ 
trast between the cores and medium is ~ 10. The influx 
of the UV photons from an emerging massive star ion¬ 
izes the tenuous warm medium within a pc from the star 
and raises the medium’s temperature to Text ~ 9000 K. 
In contrast, the temperature within the cores remains at 
Tint ~ 20 K (see Fig[T]). 

The ionization front propagates through the warm 
medium more rapidly than the sound speed. Conse¬ 
quently, the increase in Text leads to an increase in Text 
by a similar factor before the medium’s density can read¬ 
just to a new pressure equilibrium. Due to the combined 
effect of the temperature and pressure increase in both 
the warm medium and the cold medium, the Bonnor- 
Ebert mass d ecreases by a fac tor of 2.38, shifting from 
around 2 Mq (jLada et al.ll2008ll to 0.84 Mq in our model. 

^ Note that a different mean molecular weight /r is used for the 
two cases. 


Due to this sudden increase of the external pressure, 
global, coordinated star formation is induced, leading to 
a star burst. Several factors may introduce some disper¬ 
sions in the value of the modified ttibe- For example, 
the flux of ionizing photons emitted by the first 0/B 
star would be reduced by an order of magnitude if its 
mass is halved (the one we have considered in Figure [I] 
is 40 Mq). Nevertheless, the final temperature of the 
ionized region and the core only change slightly. The 
Bonnor-Ebert mass associated with slightly lower Text 
(8500 K) and Tint (15 K) would reduce mBE to 0.47 Mq. 
In this case, the Strdmgren sphere around the star will 
be much smaller, which would slow down the global star 
formation process and extend the age spread in this neb¬ 
ula. If the global star formation timescale is comparable 
or longer than the timescale of core coagulation, the dy¬ 
namics of cores may affect the IME further. Here, we 
neglect this effect on the age spread and the effect from 
the evolution of cores. Once a core becomes gravitation¬ 
ally unstable, it collapses on a free fall timescal^l. 
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(5) 

This timescale is much shorter than that associated with 
the dynamical evolution of the dense cores prior to the 
formation of the first massive stars, which is about 5 Myr 
from H13. After the external medium is ionized by the 
first massive stars, the compression of the cores reduces 
the cross section and enhances the density contrast be¬ 
tween the core and the external medium.This increases 
both the coagulation and the fragmentation timescales. 
Therefore, the further dynamical evolution of the CMF 
can be neglected. 


3. EROM THE CME TO THE IMF 

The initial stellar mass function is determined by the 
induced collapse of the dense cores with a preset CMF. 
The calculations presented here follow our previous work 
on the CMF (H13). In the previous paper, we discussed 
the evolution of the dense core mass distribution with a 
modified coagulation equation. The CMF acquires a sta¬ 
ble shape which resembles a typical observed CMF after 
several million years of evolution time. During this stage 
prior to the triggered star formation we discuss here, 
a typical star formation timescale is around 100 Myr, 
whic h keeps the starless nature of the system. Figure 
2(a) shows the pre-stellar core mass di stribution com- 
pared with the Chabrie r System IMF (iChabried l2003fl 
and recent observations (iRathborne et al.ll2009f) . The in¬ 
termediate mass range of the CMF (0.3-0.5 Mq) can be 
also parameterized with similar two power-law slopes as 
the Kroupa IMF, only shifted to a h igher mass range by a 
factor of about 3 (|Lada et al.ll2008ll . We note that both 
the observed CMF and modeled CMF have a slightly 
steeper slope than the IMF at the high mass end. How¬ 
ever, the uncertainty is also higher in those bins due to 


^ Note that after the collapse the accretion flow onto the cores 
could be still active and the star formation process may continue. 
Here, we are mainly concerned with the epoch before stellar feed¬ 
back is activated. The star formation timescale should still be 
well approximated with the free fall time tg. 
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Fig. 1.— Resulting temperature profiles in the surrounding of a 40 Mq star for different gas profiles, (a) The red line displays the 
temperature profile in a medium with constant density no = 774 cm“^. The transition from ionized to un-ionized medium happens at 
0.7pc (the Stromgen radius), (b) This setup represents a cold core at a distance of 0.6 pc from the central star. The initial density profile is 
set to no = 774 cm“^ (inter-core medium) inside 0.5 pc and no = 7300 cm~^ (core medium) beyond 0.5 pc. The solid vertical line indicates 
the inner side boundary of the core, i.e. the location of the density jump. The dashed vertical line indicates the extent of the core from 
the solid vertical line. For computational simplicity, we do not calculate the region behind the core, as the focus of this work lies on the 
penetration depth at the front side. 


small number statistics. We now take the modeled CMF 
(H13) as the initial condition and assume the global star 
formation is triggered simultaneously. During a time step 
St, the number of stars in the mass range m m + dm 
increases by 


An*{m*) 

At 



mcn{mc)r]snsf 


P{mc,m*) 

m* 


dmc 


( 6 ) 


tion of core mass finally remained in the stars, for one 
particular core. While the more often loosely used phrase 
‘star formation efficiency’ should refer to the global ratio 
of total stellar masses over total progenitor cor e mass es, 
in a certain region and time-scale. In Figure |2(b)| we 
show the IMFs generated with this model. We assume 
only cores with mass exceeding the Bonnor-Ebert mass 
can form stars, such that the star formation rate 


where n{mc) is the number density of cores within the 
mass interval {me,me + dme), f?sf(wc) is the retention 
factor (see section 3.1), 7 sf(mc) is the star formation rate 
for a dense core with mass TOc, and P{me,m*) is the 
percentage of the core mass (me) transferred into this 
specific star mass (m*) bin. 

The real star formation process may depend on differ¬ 
ent initial conditions of the progenitor cores. Here, we 
present several limiting cases. To quantify our results, 
we fit our computed IMFs with a Chabrier-like function 
(j){fj,,a,mo,^), where mo is the transition mass separat¬ 
ing the two parts of the piece-wise Chabrier-like function. 
For masses smaller than mo, the model IMF is described 
by a lognormal form with mean around ^ and disper¬ 
sion a; while the distribution at masses larger than mo 
is described by a power-law with index 7 . We require 
the distribution to be continuous at mg. The best fitting 
parameters are shown in Tabl e [1] We compy e these fit¬ 
tings with the Chabrier IMF (lChabrieill2003ll model and 
discuss the limits of different prescriptions. 

3.1. Case 1: burst of single stars 

Given the similarity between the CMF and the IMF, 
some authors suggested an one -to-one conversio n of 
dense cores into young stars fe.g.. lLada et al.l[2008[ l. Al¬ 
though this scenario has been cri ticized on the basi s of 
binary stars’ prevalence (see e.g.. iSmith et al.ll2009ll . it 
is nonetheless informative to explore this simplest possi¬ 
bility. In this case, P{me, m*) = 6{m* — mc? 7 sf (me)) (see 
Equation| 6 |). rjsf is the retention factor, which is the frac¬ 


7sf(™c) 


Jo if me < msE 
\Fsf if me > mBE 


(7) 


Here, Fgf = tf^^ (see Equation [S]) is the characteristic 
star formation rate. We also assume that the retention 
factor, i.e. the amount of the core mass retained in the 
resulting star, is assumed to be constant with a value 

0 < 77sf(mc) < 1. 

The results in Figure |2(b)| show that the IMFs are 
shifted toward the low-mass end, with a broad peak 
near the original Bonnor-Ebert mass modified by the 
retention factor. With rjsi = 0.3 the new peak is at 
around 0.7 Mm, w hich corresponds to Lada’s suggestion 
(|Lada et al.ll2008f) . The overall shape of the IMF resem¬ 
bles the observed stellar IMF. The low-mass end have 
similar dispersion with that of Chabrier’s IMF, provided 
the retention factor r]sf is not a sensitive function of the 
progenitor core mass. The high- mass end powe r-law is 
steeper than the Salpeter slope ()SalDeteilll955fl . where 
the observational uncertainty in this range is also rela¬ 
tively large. However, the minimum mass of stars we can 
produce with a single star formation prescript, is depen¬ 
dent on the critical Bonnor-Ebert mass, and rjsf such that 
Wcut-off ~ WBE,new 7 sf- Due to this cut-off cvcn with a 
low value of r/sf (r/sf =0.3, dark grey line in Figure|2(b)|), 
this model has no inference for the lowest mass stars. 
The bulk of the modeled IMFs are more massive than 
the Chabrier IMF. We conclude that a non-unity reten¬ 
tion factor alone would not explain the transition from 
CMF to IMF. 
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Fig. 2.— CMFs and IMFs generated in our model under different assumptions, (a) The CMF obtained from H13 (solid black li ne) is 
adopted as the initial condition for our calcula tions. Grey line shows the observed CMF in the Pipe from iRathborne et al.l II2009I) . For 
comparison, the Chabrier IMF (|Chabriei1l2Q03l') . shifted to the low mass end by a factor of 3, is shown as the dashed blue line, (b) Grey 
curves are displaying the stellar IMFs for single stars with different retention factors T/gf (me). Any core evolution through coagulation or 
fragmentation is neglected, (c) The black solid line is showing the generated total IMF if all the cores undergo binary fragmentation with 
(3 = 0.2 (see Equation 0 and 77sf(mc) = 0.3. (d)The black solid line is showing the total IMF, assuming 100% binary star formation. 
The secondary stellar masses follow a power-law dis tribution with index q = 1.5 (see Equation [Oj . We show both the Chabrier IMF (blue 
dashed curve) and the Kroupa IMF JKroupal 120021) (grey dash-dotted lines) as a comparison. The modeled IMF is fitted with a broken 
power-law (red dash-dotted lines). 


TABLE 1 

IMF FROM DIFFERENT STAR FORMATION MODELS. 


distribution 

M{Mq) 

mean 

^(Mq) 

dispersion 

mo{MQ) 
transition mass 

7 

power-law index 

SSR 

Sum of Squared Residuals 

Chabrier 

0.22 

0.57 

1.00 

1.30 

N/A 

singleO.3 

0.66 

0.54 

4.05 

2.19 

0.000011 

singleO.5 

1.05 

0.54 

6.42 

2.18 

0.000011 

singleO.7 

1.55 

0.54 

9.36 

2.14 

0.000011 

singleO.9 

1.95 

0.55 

11.02 

1.98 

0.0000097 


Note. — List of parameters for the resulting IMFs modeled with a Chabrier type function. We 
fitted the resulting IMF in the mass range 0.1-30 Mq with m ^ ^ ^ with 

dlog m foL > ™0- Continuity is imposed at mo. Columns starting with ‘single’ refer to the 

results from the single star formation model in section 3.1, and the number attached to it corresponds 
to the retention efficiency. 
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Fig. 3.— Analytical results for f(q), i.e., the probability distri¬ 
bution function (PDF) of observed mass ratios {q = msec/niprim) 
in binary systems. The conditional distribution of mass ratio given 
the progenitor core mass is assumed t o b e a uniform distribution 
in regardless of mcore (see Ea ll3ll . For the CMF prescrip¬ 

tion dn/dm = m~°‘, the blue line corresponds to o = 1.3 and the 
red line corresponds to a = 2.35. The total accumulated prob¬ 
ability is renormalized to 1 and = 0.05. The dashed line 

is the best fit power-law ( f(q) oc to the observations from 

IKouwenhoven et al.l II2007I1 . 


3.2. Case 2: Solely binary formation 

A large fraction of the young stars are binaries. Each 
component of the binary contributes to the statistics of 
the IMF. Another limiting case of interest is the possibil¬ 
ity that all the cores with mass in excess of the modified 
Bonnor-Ebert mass would fragment into binary stars. In 
principle, binary star formation is a consequence of ro¬ 
tational fragmentation and the kinematic properties are 
determined by the angular momentum distribution of the 
cores and their cooling ability. However, we are primar¬ 
ily interested in the IMF rather than the period distribu¬ 
tion. We adopt an idealized power-law distribution for 
the mass ratio q such that 


fin* 

“'‘binary 

dq 


OC q 


( 8 ) 


Here, the cut-off is taken as = 0.05. We explore a 
range of power-law indices j3 between 0.2 and 2.5, which 
essentially covers the extreme limits. We assume a con¬ 
stant retention rjsi = 0.3 and star formation rate same as 
in Equation 0 . The results here remain valid if the re¬ 
tention factor does not vary strongly with the progenitor 
core mass. 

All the resulting IMFs with different /3 display a broad 
shifted peak arou nd the new, reduced Bonnor-Ebert 
mass (Figure 2(c)). The results are very slightly affected 
by different pm this wide range, so only P = 0.2 is 
shown. The mass associated with the turn-over of the 
IMF is smaller than that found for the single star forma¬ 
tion model. The segments in 0.46 Mq < m < 1.45 Mq 
with power-law index 71 = 0.489 and in 1.45 Mq < m < 
10 Mq with power-law index 72 = 1.5. The two indices 
are quite similar to those of the piece-wise Kroupa IMF. 
We expect that the allowance of trinaries and quaternar¬ 
ies will further move the peak towards low-mass end. 


In a classical study of the bin ary star population census 
bv iDuauennov fc Mavod (11991! ) . the IMF for the compan¬ 
ions of G-dwarf stars is thoroughly analyzed. They show 
that binaries on average can be formed by random combi¬ 
nation of stars drawn from the same IMF. Following their 
basic approach, we consider the IMF of the primary and 
secondary stars separately. We assume all the progeni¬ 
tor cores fragment into binaries and the secondary star 
masses follow a distribution (e.g. Gaussian. Salpeter or 
Miller-Scalo (jMiller fc Scalolll979D power-law) indepen¬ 
dent of their primaries’ mass. For simplicity, we adopt a 
power-law, so the transfer functions as in Equation ([ 6 ]) 
become: 


.^sec(^^sec; ^^c) OC dnsec/drUgec OC (ms0(;) , (9) 


.Pprim(^sec; ^c) — Ps&cijl^c ^secj^c); (1^) 

where the subscript ’sec’ and ’prim’ represent the sec¬ 
ondary and the primary star. For Pggc the range of 
TO*is 0.08 Mq < msec < 0.5 rjsimc. The power-law in¬ 
dex for the secondary star mass distribution, a, varies 
from 0.8 to 1.5 in our calculations, while the star forma¬ 
tion rate is as in Equation 0 and the retention factor 
rjsi = 0.3 is still assumed. Here, the IMF of the primary 
star preserves its characteristic broad peak and overall 
shape. The characteristic mass associated with this peak 
is smaller than that of the CMF. 

The consequent IMF for the primary stars is very 
slightly influenced by a in this wide range, so w e only 
show the result with a = 1.5 in Figure |2(d)| Fol¬ 
lowing [Duquenno^&iS^^ ()1991l) . we compare the to¬ 
tal (primary and secondary) simulated IMF with some 
well-known IMFs, such as the Kroupa IMF. The shape 
of the modeled IMF resembles the Kroupa IMF more 
closely at masses below 0.7 Mq. We fit it with two 
broken power-laws. The low-mass end (m < 1 Mq) 
completely reproduces the assumed power-law we put in, 
with dn/dlog m = while the intermediate mass 

range (1 Mq < m < 10 Mq) recover the Salpeter power- 
law dn/dlog m = 

3.4. Analytical results for the binary mass ratio 
distribution 

An alternative approach to characterize binary star 
statistics is to utilize the mass ratio, defined as q = 
rnsec/’Tiprim- Observations show that the mass ratio also 
roughly follows a power-law distrib ution, as in Scorpius 
OB2 for intermediate mass stars (jKouwenhoven et al.l 
120071) . Given a probability distribution of mass ratio 
p(g|mcore) in a binary formation process (formed from 
cores with mass mcore), we can combine it with the CMF 
piiTicore) to predict the mass ratio distribution f{q): 

^TTlmax 

p{q) = / p(g| mcore )p (mcore )dmcore- (H) 

d mmin 

With a power-law mass distribution of the progenitor 
dense cores and an assumed constant retention, we can 
estimate the probability of finding a binary system with 
stellar masses of mprim and msec to be 


3.3. Case 3: The binary companion’s IMF 


/(mprim,m sec) CX ('m-priin “^sec) 


( 12 ) 




























7 


where a is the power-law index of the progeni- 
tor CMF. For random p airing of binary systems 
(|Duauennov fc Mavoi]ll99lli . 

Pq\mcoTe) = -r-^ - iorq G 1] (13) 

^ ^min 


with Qjjjjjj as the minimum allowed binary ratio. 

The cumulative probability P{q < x) for all the binary 
systems would then be 


rXTYlrn.a.x 

P{q < x) = i dm^ 


dlTlsecf (^prim; ^sec)- 


(14) 

Thus, the statistically averaged distribution of mass ra¬ 
tios q would be 


tt-,-dP{q) l + {2-a)q+{l-a) 

J{q) = —;— oc 


»min i 


dq 


(2 - a)(l -b qY 


-(!+<?)' 
(15) 

where mniax(= Vsimcore) is the upper limit for progen¬ 
itor core mass. We choose a values corresponding to 
Kroupa’s IMF power-law indices. Although there are 
some uncertainties in the minimum allowed mass ratio, 
binary systems with q ratio around 0.05 or less has been 
observed so we adopt = 0.05 in the calculation. 
Results are shown in Fi gure B] and compared wit h the 
observed power-law from lKouwenhoven et al.l (1200711 . Al¬ 
though the distribution of q depends on the value of Qjjjjn, 
the power-law slope converges, for q» 0.05, to the ob¬ 
served value. 


4. DISCUSSION AND CONCLUSIONS 

In this work, we continue our investigation on the ori¬ 
gin of the stellar IMF. Based on the similarity between 
the observed stellar IMF and the CMF of dense star¬ 
less cores in the Pip e nebula, we assume that they are 
closely connected. In iGritschneder fc LinI (1201211 . we sug¬ 
gest that the cold cores of molecular gas are the byprod¬ 
ucts of a thermal instability or the fragmentation of the 
a shocked shell and they are pressure confined by ten¬ 
uous warm atomic medium. The CMF of these cores 
is the natural outcome of their collisional coagulation 
and their fragmentation due to their hydrodynamic in¬ 
teraction with the surrounding medium (HI 3). We also 
assume that these cores become unstable, undergo grav¬ 
itational collapse, and evolve into protostars after their 
mass exceeds the Bonnor-Ebert mass. 

Although this simple model reproduces the basic ob¬ 
served slopes of the stellar IMF, there is a factor-of-three 
offset between the stellar mass associated with the peak 
of the stellar IMF and that associated with the peak of 
the cores’ CMF. The main motivation of the investiga¬ 
tion presented in this paper is to suggest a mechanism 
to account for this shift. 


In iHuang et al.l (|2013[1 . we suggest that the peak of the 
CMF is essentially the Bonnor-Ebert mass of the cores. 
In typical molecular clouds such as the Pipe nebula, the 
Bonnor-Ebert mass is > 1 Mq. In principle, cores less 
massive than the Bonnor-Ebert mass are stable and they 
do not turn into stars. Yet, many low-mass stars are 
formed in young stellar clusters. 

In this work, we make an attempt to resolve the issues 
of 1) the offset between the mass associated with the peak 
of the CMF and the IMF, and 2) the prolific production 
of sub-solar type stars in molecular clouds. We demon¬ 
strate that the onset of first massive stars in these clouds 
photoionize and heat the surrounding medium without 
significantly changing the ionization fraction and tem¬ 
perature of the cores. This feedback effect largely in¬ 
creases the external pressure which confines the cores 
(by up to two orders of magnitude). This change leads 
to the compression of the cores and a reduction in their 
Bonnor-Ebert mass. 

The collapse of cores with a mass greater than the 
modified Bonnor-Ebert mass (and less than the origi¬ 
nal Bonnor-Ebert mass) can now lead to the formation 
of a large population of sub-solar-mass stars. The stellar 
IMF generally preserves the shape of the Salpeter-like 
CMF with a significant lowering in the peak mass. The 
shape of the IMF (dispersion and high mass end slope) 
is not strongly modified by the formation of binary stars 
through rotational fragmentation. In our model, the in¬ 
clusion of binary fragmentation during the collapse is es¬ 
sential for the production of stars with mass lower than 
the modified Bonnor-Ebert mass. 

One implication of this induced star formation scenario 
is that the intrinsic age spread of the stellar cluster is 
naturally very small. In our analysis, we have adopted 
an idealized treatment of the retention efficiency (often 
loosely referred to as star formation efficiency). A mass 
dependence in the retention factor may lead to some 
modification in the IMF. We present several models for 
single and binary star populations. They generally repro¬ 
duce the transition from the CMF to the IMF suggested 
by the observational data. A prolific production of triple 
and hierarchical systems may also lead to the formation 
of very low-mass stars. 

Finally, these models are constructed with solar com¬ 
position. In metal-deficient gas, such as protoglobular 
cluster clouds, the inability to cool may lead to a higher 
internal temperature and a higher Bonnor-Ebert mass 
in the cores. This feedback mechanism may be particu¬ 
larly important in triggering the formation of low-mass 
stars with lifespans in excess of 10 Gyr. We will further 
explore this possibility in the future. 
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